Directed motion and useful work from an isotropic nonequilibrium distribution 
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We demonstrate that a gas of classical particles trapped in an external asymmetric potential 
undergoes a quasiperiodic motion, if the temperature of its initial velocity distribution Tm differs 
from the equilibrium temperature, Teq. The magnitude of the effect is determined by the value 
of Tne — Teq, and the direction of the motion is determined by the sign of this expression. The 
"loading" and "unloading" of the gas particles change directions of their motion, thereby creating a 
possibility of shuttle-like motion. The system works as a Carnot engine where the heat flow between 
kinetic and potential parts of the nonequilibrium distribution produces the useful work. 

PACS numbers: 05.20.-y, 05.40.-a, 05.40.Jc 



Since the Maxwell demon thought experiment, the ex- 
traction of useful work and directed motion from unbi- 
ased nonequilibrium distributions has been the source of 
fascination, intrigue and confusionji Being a fundamen- 
tal scientific problem, it is also of significant practical 
interest for various biological and nanotechnological ap- 
plications. Biological systems driven by the nonequilib- 
rium fluctuations include translocases, which pull pro- 
tein across membranes; kinesin, which walks along micro- 
tubules; myosin, which moves along actin in muscle; heli- 
cases, which unwinds DNAfS Recently, there have been a 
flurry of activities to harness nonequilibrium fluctuations 
to create various types of working nanodevicesi^i^ In this 
Letter, we propose a new type of "motors" which are 
driven by the heat flow between nonequilibrium velocity 
and equilibrium coordinate distributions. 

Consider a gas of classical non-interacting particles 
trapped in external potential U{x). If the gas is in ther- 
mal equilibrium with a heat bath at the temperature T^q, 
then positions {x) and velocities (v) of the gas particles 
have the equilibrium Boltzmann distribution: 

Pb{x,v) ^ Pb{x)pb{v); (1) 

Pb{x) — exp{ — [/(x)}/ j c?x exp{— J7(a;)}, 

Pb{v) = v/l/(27r)exp{-«V2}. (2) 

We use the reduced variables in this Letter: Time, po- 
sition, velocity, and energy are measured in units of 
y/mL'^/{kBTeq), L, ^ykBTeq/m and kBT^q, respectively 
{kB is the Boltzmann constant, m is the mass of a gas 
particle and L is a unit of length) . 

Since equilibrium fluctuations cannot produce any di- 
rected motion, the average velocity of the gas {v{t)) g — 
J dvdxpB(v, x)v{t) is identically zero at any time t. What 
happens if we instantaneously change the velocity distri- 
bution of the particles to some non-equilibrium distri- 
bution Pne{v)l Apparently, if Pne{v) — 6{v — v'), then 
the gas undergoes a periodic motion within the potential 
well. Such a distribution (and akin non-singular distri- 
butions) are biased, since they predict nonzero mean ini- 
tial velocity {v{t = 0)) = v' . In this Letter, we consider 



more physically interesting and experimentally accessible 
situation, when the nonequilibrium distribution is unbi- 
ased and depends on the modulus of the particle velocity, 
Pne{\v\), SO that (t)^"+^(t = 0)) = for any integer n. To 
be more specific, we assume for simplicity that Pne(|t'|) 
can be represented by the Boltzmann distribution, but 
at a nonequilibrium temperature T„e ^ Tg^: 



Pne{x,v) = pB{x)pne{v), 



1 T, 



Pneiv) = \/ 7^77^ exp <^ - — — 



27rT, 



2 T„ 



(3) 



Such nonequilibrium distributions over molecular ve- 
locities and angular momenta are produced, for exam- 
ple, in the course of photodissociation A + Tiv ^ B + 
products,— i^"^!^ which is normally considered as instan- 
taneous on the timescale of molecular translation and 
rotation. "'^'^ Indeed, let x^, p,;, and be the positions, 
the momenta, and the masses of the parent {i — A) and 
product (i — B) molecules. Since the parent molecule 
falls apart instantaneously, the positions do not change, 
= x^. As to the momenta, we can write^ = 
(mB/™A)Pyi + ^- The first term here describes mapping 
of the parent molecule translation into that of the prod- 
uct. The quantity A depends on the carrier frequency 
V of the dissociating laser pulse and originates from the 
impulsive force arising due to the rupture of chemical 
bond(s) of the parent molecule. If the parent molecules 
possess initial equilibrium distribution over their posi- 
tions and momenta, then dissociation preserves the orig- 
inal distribution over the positions, /9b(x), but produces 
a nonequilibrium distribution over the momenta of pho- 
toproducts, p„e(p)- If A 7^ 0, then p„e(p) is a shifted 
Gaussian distribution. If the laser is tuned as to make 
A = 0, then p„e(p) is the Boltzmann distribution at 
the nonequilibrium temperature T„e = (mB/rnA)Teq. In 
general, any reaction of the kind A ^ A' (photoassoci- 
ation, photoisomerization) will create a nonequilibrium 
distribution over velocities of the species A and A' , pro- 
vided the processes is rapid enough on the timescale of 
molecular dynamics. 

We begin with the standard Newton equations of mo- 
tion for non-interacting particles in potential U{x): 



V = X, i) = —dU{x)/dx. 



(4) 



In the equilibrium case, initial values of positions and 
velocities of the gas particles should be sampled from the 
Boltzmann distribution ([1]). Since our aim is to evaluate 
the average velocity, 



v{t)) = / dvdxpne{v,x)v{t), 



(5) 



we need to solve Eq. ([4]) with the initial nonequilibrium 
distribution ([3|). 

This average velocity {v{t)) possesses several rather un- 
usual properties. It is well known that any equilibrium 
correlation function calculated for any Hamiltonian sys- 
tem is an even function of time. This is a direct conse- 
quence of the time- inversion symmetry. In contrast, as is 
easy to demonstrate, 



2n+l 



n=0 



{2n + l)\ 



(6) 



is an odd function of time, so that {v{t)) = — {v{—t)). 
This type of behavior is caused by the special initial 
nonequilibrium preparation of the gas at i = 0. We can 
calculate the first few terms in series ([6]) analytically: 



Go — 0, Gi 



T 



q \ dx^ 



(7) 



As one may expect, Gi is proportional to the difference 
between the equilibrium and nonequilibrium tempera- 
ture, since all G„ must vanish at equilibrium. Eq. ([7]) 
shows that the external potential must not have a sym- 
metry axis, in order to get {v{t))^ 0. 

To be more specific, we consider a gas trapped in the 
nonlinear asymmetric potential 



U{x) = ^^2,2/2 + a{x ~ s)3 -t- b{x - s)-^ 



(8) 



Here lo is the harmonic oscillator frequency and the pa- 
rameters a, 6, and s describe the nonlinear part of the 
potential. One cannot analytically calculate {v{t)) for 
this potential. Therefore, we will compute it numerically 
but, first, it is instructive to show some results of the per- 
turbation theory. Assuming that the anharmonic part of 
the potential is small and retaining linear in a and b con- 
tributions only, we get 



T —T 

J- ne e 



eq 



(3a- 12s5)$(i), 



eq 



where 



2 



sin((jjt){l — cos(ti;t)}. 



(9) 



(10) 



So, {v(i)) exhibits a periodic motion. Its magnitude 
and direction are determined by the value and sign of 
Tne — Tf,q. Sincc the difference in temperatures is readily 
linked to the difference of masses in distribution Pne{v), 
we see that "loading" and "unloading" of the particles 
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FIG. 1: Mean velocity {v(t)) for ideal gas of particles moving 
in the potential (|8]) with u) = 1, s = 0.2, a = 0.003, and 
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FIG. 2: Same as in Fig. 2 but for s = 1, a = 0.3, and h = 0.5. 



(that is the increase and decrease of their masses) changes 
directions of their motion, thereby creating a possibil- 
ity of shuttle-like motion. Furthermore, if we consider 
dissociation A ^ B + C , then (different) particles B 
and G are produced at different nonequilibrium temper- 
atures and thus induce different fluxes, allowing for the 
mass separation. ^-^ The terms stemming from the cubic 
and quartic branches of U{x) enter Eq. ^ with dif- 
ferent signs and thus lead to the opposite senses of mo- 
tion. Interestingly, the contributions cancel each other 
for a = Ash, but this is valid only within the perturba- 
tion theory. Higher order terms ensure nonzero {v{t)). 

We also compute {v{t)) exactly, solving numerically 
equations of motions (jH) by Verlet algorithm for differ- 
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FIG. 3: Mean velocity {v{t)) for ideal gas of particles moving 
in the potential ([8]) with uj = 1, s = 5, a = 0, and 6 = 5; 



FIG. 4: Total ensemble-averaged work {v{t)^ — '^(O)'^) for gas 
of particles moving in the potential ((Sjl with cu = 1, s — 1, a — 
0.3, and b = 0.5; for Tn^ = 2Teg (full line) and r„e = Teg/2 



ent initial conditions sampled according to the noncqui- 
librium distribution ([3]). Typical behaviors of {v{t)) are 
presented in Figs. 1-3. If nonlinear corrections to U{x) 
are small, then several first oscillations of {v{t)) in Fig. 
1 are excellently described by the perturbative formulas 
(HI) and ([TU)) . However, the oscillations decay in time, de- 
spite absence of dissipation. This can be understood by 
switching to the action (/) - angle {ip) variables. Then 
we can write 



/ dIdLppne{I,f)vn{I)exp{in{if~uj{I)t)}, 



(11) 

Vn{I) = (27r) ^ J d(pv{I , Lp) e-xjp{inip} (evidently, vo{I) = 
0). As distinct from the harmonic oscillator, the fre- 
quency to{I) of a nonlinear oscillator depends upon the 
action variable and {v{t)) decays due to dephasing. 

If the the nonlinear contributions to U (x) become more 
pronounced, then the magnitude of {v{t)) increases (Figs. 
2 and 3), but the oscillations decay more rapidly. It is 
remarkable that {v{t)) can exceed the mean equilibrium 
thermal velocity ^JkBTeq/m (Fig. 3). The corresponding 
value of the nonequilibrium temperature T„e = SOTeq is 
not unreasonable, e.g., for hot CN fragments produced 
through photodissociation of ICN (see Ref<^). 

The quantity (| j^^^ F{x)dx^ = {v{t2f - v{tif) /2 can 

be regarded as the total ensemble-averaged work per- 
formed by the gas. Evidently, this quantity is identically 
zero at equilibrium. Under nonequilibrium conditions, 
this is nonzero oscillatory function (Fig. 4) . The presence 
of oscillations means that, by choosing the appropriate 
time moments ti and t2 , we can make the averaged work 
{v{t2y — w (ii)^) /2 positive, irrespective of whether the 
nonequilibrium temperature T„e is smaller or higher than 
the bath temperature T^q. So, our system can be re- 
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FIG. 5: Mean velocity {v{t)} for a moderately damped (C = 1) 
Brownian particle moving in the potential ((Sjl with u = 1, s — 
1, a — 0.3, and b = 0.5. Full line corresponds to r„e ~ 2Teq 
and dashed line corresponds to r„e — Teq/2. 



garded as a true heat engine, which produces a cycled 
motion out of unbiased nonequilibrium distributions and 
can be harnessed to make a useful work. 

Now we turn our attention to the role of the system- 
bath interaction. We thus convert the Newton equation 
of motion Q into the Langevin equation: 



V = X, V 



-dU{x)/dx-iv + f{t). 



(12) 



Here, ^ is the friction and f{t) is the stochastic white- 
noise force, which obeys the fluctuation-dissipation rela- 
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tion {fit) fit')) = 2^S{t - t'). The coupling to the heat 
bath restores equihbrium distribution ([1]), so that eventu- 
ally {v{t)) when t — > oo. As in the bath-free case we 
can obtain the perturbative solution of Eq. p^ . which 
is valid up to the terms linear in a and b. The result is 
again given by Eq. ([9|), in which one has to put 



2 , 



— cx-p] — — r I 



D 



exp{-CO 



exp{VDt} 



exp{-\/iJ<} 



(13) 



D = — Auj^. In the underdamped case {D < 0), 
<i>(t) consists of two contributions, which oscillate at fre- 
quencies \J—D and y/—D/2, and decay ~ exp{— ^t} and 
exp{— correspondingly. Furthermore, there ex- 
ists a monotonically decaying contribution ^ exp{— ^i}, 
which is absent in the dissipation-free case. If 2uj ^ ^, 
then Eq. reduces to the damped analogue of Eq. 

m 



'^{t) = ^sin(wi)exp{--t} ( 1 - exp{--t} cos(wt) 



(14) 

We solve Langevin equation (|12p numerically, as de- 
scribed iuji^ Fig. 5 depicts {v{t)) calculated for mod- 
erate damping ^ = 1. A comparison of Figs. 2 and 5 
reveals that damping, as expected, decreases the abso- 
lute value of {v{t)), but the effect definitely persists. The 
damping does not destroy the oscillation of the ensemble- 
averaged work shown in Fig. 4. This means that the phe- 
nomenon of {v{t)) 7^ is quite robust and may exist in 
dense gases and even liquids. In the overdamped case 

> 2a;, i > 1/^) $(t) becomes - 



m = ^exp{-— t} (^2exp{-— 



1 



(15) 



Note, however, that the question about the existence of 
{v{t)) 7^ is not only about the strength of friction, but 
also about the corresponding timescale. Since motion of 
any Hamiltonian system is dissipation-free for t 1/^, 
{v{t)) 7^ on that time interval. 

It has not escaped our notice that, with suitably se- 
lected external periodic potential, the system represents 
a new type of motor which convert the nonequilibrium 
disparity in the velocity and coordinate distribution into 
directed motion. Although both standard Brownian 
ratchets^^'A^t!^*' and our motor use the broken symme- 
try of the external potential and nonequilibrium fluctu- 
ations to produce the directed motion, there are three 
important differences between them. First, the standard 
Brownian motors require external time-dependent (de- 
terministic and/or stochastic) forces lacking the detailed 
balance symmetry and perpetually driving the system 



out of equilibrium. Alternatively, the temperature can 
be either driven externally or kept position-dependent. 
In our case, no external driving or nonuniform temper- 
ature are necessary. The directed motion is entirely de- 
termined by the nonequilibrium system preparation, and 
it is not essential if the initial (nonequilibrium) kinetic 
energy fcBT„e/2 is greater or smaller than the equilib- 
rium kinetic energy, ksTeq/i. Second, with a few notable 
exceptionsii^^i^OeiiOsag^^^^^^^ the Brown- 

ian motors are normally studied in the overdamped limit. 
In any case, the inertia is not a key factor ensuring the 
directed motion of Brownian motors (although it is essen- 
tial for the phenomenon of absolute negative mobilityi^ 
and sometimes is responsible for the very occurrence of 
the directed transport^^). On the contrary, our motor is 
operating in the inertial regime, and the effect is more 
pronounced in the gas phase, when friction ^ ^ 1 (Eq. 
(fT4|) ). The effect does not disappear in the overdamped 
limit, but becomes very small, {v{t)) ~ (Eq. (fT5|) ). 
This is entirely understandable, since the higher is the 
friction, the more rapidly Pne(v) relaxes to pb(v). How- 
ever, if a reversible reaction A < — > A' takes place in 
our ensemble, then Pne{v) can be recovered periodically. 
Third, the Brownian motors can be envisaged as been 
connected to two different heat bathes.— i^^i^i Similarly, 
the "blowtorch effect" (i.e. a nonuniform temperature 
profile along the potential surface, which, in the simplest 
case, can be realized within a two-bath model) can also 
induce a directed current i^i^S In our case, we can think 
of our motor as being connected to kinetic and potential 
bathes, which are initially at different temperatures. In 
this sense, it can be treated as a molecular cycled Carnot 
engine. 

Note that nonlinear dissipation effects can produce a 
directed motion of the Brownian particle provided its ini- 
tial velocity distribution is asymmetric. '^^ 

The results presented in this work are not limited to a 
specific Gaussian form of the nonequilibrium distribution 
Any distribution p„e(|t;|) which differs from pb{v) 
gives rise to nonzero and periodic {v{t)). Realistic veloc- 
ity distributions for dissociating molecules are described 
in papers i^'^'^" The theory can also be straightforwardly 
extended to rotational and vibrational motion. 

Summarizing, we demonstrate that the gas of non- 
iteracting particles trapped in an external asymmetric 
potential moves with nonzero oscillatory mean velocity 
{v{t)), provided its initial velocity distribution is pro- 
duced at a temperature r„e which differs from the equi- 
librium bath temperature, Tgg. The magnitude of the 
effect is determined by the value of T„e — Teg, and the 
direction of motion is determined by the sign of this ex- 
pression. Practically, such a nonequilibrium distributions 
over molecular velocities and angular momenta are pro- 
duced in the course of photodissociationf^ii^i^ or due 
to any reaction of the kind A A' , provided it is rapid 
enough on the timescale of particle dynamics. In molecu- 
lar dynamics simulations, we can instantaneously rescale 
molecular masses (that is equivalent to the rescaling of 
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temperature) and arrive at the nonequilibrium distribu- 
tion ([3]). The phenomenon of occurrence of {v{t)) ^ 
is robust and does not require any careful tuning of the 
of the underlying parameters. It is not restricted to one- 
dimcnsional case and any asymmetric potential produces 
the effect, which will remain in the quantum case, too. 
Furthermore, our considerations do not contain any spe- 



cific timescale and lengthscale, and are valid both for 
atomistic ensembles and for collections of macroscopic 
objects. 
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